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Abstract. We note that the solution of accretion disk dy- 
namics for an initial delta-function mass distribution gives 
a light curve that fits both the rise and the decay pat- 
tern of the outburst light curves of black-hole soft X-ray 
transients (BSXTs) until the onset of the first mini out- 
burst quite well. The Green's function solution of Lynden- 
Bcll & Pringle (1974) is employed for two different time- 
independent viscosity laws to calculate the expected count 
rates of X-ray photons in the Ginga energy bands as a 
function of time. For both models basic characteristics 
of the outburst light curves of two typical sources GS 
2000+25 and GS/GRS 1124-68 are reproduced together 
with plausible values of the thin disk parameter a and 
the recurrence times. This agreement with the outburst 
light curves and the source properties during quiescence 
support the idea of mass accumulation and the sporadic 
release of accumulated mass at the outer disk. 

Key words: Soft X-ray transients, accretion discs - 
black-holes, X-rays 



1. Introduction 

About two thirds of the known LMXBs are persistent and 
the remaining one third are transient sources (van Paradijs 
1995). The transient sources exhibit a soft X-ray spectrum 
during outbursts (soft X-ray transients — SXTs). An opti- 
cal outburst also occurs together with the X-ray outburst. 
During an outburst the X-ray and the optical properties 
of black-hole soft X-ray transients (BSXTs) and neutron 
star soft X-ray transients (NSXTs) are very similar to each 
other and also to those of persistent sources. X-ray lu- 
minosities of SXTs increase from below 10 33 erg s _1 to 
~ 10 37 -10 38 erg s _1 during an outburst. The time scales 
for decline range from several tens of days to more than 
one hundred days. Although some sources show rather 



monotonic decline in many cases there is more com- 
plex behavior. Some sources remain persistent for more 
than one year. Most of the light curves show a steepen- 
ing around 10 36 erg s _1 corresponding to a mass accre- 
tion rate M x ~ 10 16 g s -1 inferred from the X-ray lu- 
minosity. Later they turn back to their quiescent states 
with L x ~ 10 31 - 10 33 erg s" 1 and M x ~ 10 11 - 10 13 
g s _1 (Tanaka & Shibazaki 1996 and references therein). 

There are nine known X-ray binaries with strong black- 
hole candidate primaries (van Paradijs 1995). Three of 
these sources are non-transient high mass systems and 
the remaining six are transient LMXBs (Table 1). Four 
BSXTs in this group, A0620-00, GS 2000+25, GS/GRS 
1124-68 and GRO J0422+32 show striking similarities. 
Their outbursts have fast rise times around a few days. 
The decay phases can be fit with exponentials with 
time constants around one month. During decline all 
four sources exhibit secondary (mini) outbursts. A0620- 
00, GS/GRS 1124-68 show also tertiary outbursts. The 
tertiary outburst is absent in the light curve of GRO 
J0422+32, but it may also be present in GS 2000+25. 
Secondary maxima were detected around 80 days af- 
ter the onset of the main outburst of GS 2000+25 and 
GS/GRS 1124-68. A0620-00 and GRO J0422+32 exhib- 
ited the secondary maxima 50 days and 125 days af- 
ter the onset respectively. For both GS 0620-00 and 
GS/GRO 1124-68 the tertiary maxima were observed 
about 200 days after the onset of the main outburst 
(Tanaka & Shibazaki 1996). The basic features of the mini 
outbursts are the sharp increase in luminosity by a factor 
of ^1 — 3 and decay patterns which mimic the decay 
after the first outburst. A0620-00 is a recurrent transient 
with a recurrence time of 58 years (Tsunemi et al 1989). 
The outbursts of the other sources were detected only once 
putting a lower limit to their recurrence times of around 
a few ten years (see Table 1). These similarities allow the 
working hypothesis that all of these sources run with a 
similar mechanism. 
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Since the conditions for outburst build up in the qui- 
escent state it is important to treat the observations made 
during quiescence. According to the standard disk model 
(Shakura & Sunyaev 1973) most of the X-rays (> 70%) 
come from the inner parts of the disk extending to a ra- 
dius which is around ten times the radius of the last stable 
orbit, R = 3R g = 3(2GM x /c 2 ). During the quiescent 
state it is possible to fit a black-body curve roughly to the 
observed X-ray spectrum. However, inner disk tempera- 
tures of ~ 0.2 — 0.3 keV obtained from these fits give very 
small values (1 — 10 km 2 ) for the X-ray emitting area of 
the disk, while realistic areas imply temperatures that are 
orders of magnitude smaller than the soft X-ray temper- 
atures that fit the spectra. A disk becomes optically thin 
below a certain accretion rate depending on the viscosity 
parameter a (Shakura & Sunyaev 1973). An optically thin 
or gray disk may be one way of explaining this inner disk 
problem. 

The mass accretion rates obtained from the optical ob- 
servations during quiescence (~ 10 15 — 10 16 g s _1 ) are or- 
ders of magnitude larger than those (~ 10 11 — 10 13 g s _1 ) 
obtained from the X-ray luminosities (see, e.g., McClin- 
tock et al 1995, for these questions in A0620-00). 

An advcction dominated inner disk has been proposed 
to account for the properties of BSXTs. According to this 
model most of the energy released in the inner disk is 
advected into the compact star instead of being radiated 
from the disk (Narayan et al 1996). Regarding the differ- 
ence between mass accretion rates obtained from optical 
and X-ray luminosities during the quiescent states, the 
NSXTs are similar to BSXTs. Although this model could 
explain the observed spectra and also the low X-ray lumi- 
nosities of BSXTs in quiescence it is not able to explain the 
properties of NSXTs because of the existence of the solid 
surfaces of the neutron stars: Whatever the disk structure 
is matter finally reaches the neutron star surface which 
should produce a distinct black-body like emission. 

Another possibility is mass accumulation in the outer 
regions of the disk. This accumulated matter may be the 
source of the outburst. For both BSXTs and NSXTs, the 
mass accretion rate inferred from optical observations may 
indicate the mass accretion rate arriving to accumulate in 
the outer disk, while the much lower M x inferred from the 
X-ray luminosity may reflect the trickle of mass that pro- 
ceeds through a mass-starved inner disk characteristic of 
the quiescent state. This qualitative picture would persist 
until the mass accumulated in the outer disk reaches the 
level critical for an outburst. 

Models for SXTs must address the characteristic time 
scales: (1) Fast rise of the light curves around a few days; 
(2) The subsequent decay with a time scale of the order of 
one month; (3) Recurrence times of the order of a few ten 
years. Further, as viscosity is involved in determining the 
decay, values of the viscosity parameter must be plausible. 
For comparison a ~ 0.1 — 1 for dwarf novae and related 



systems for the reaction time of the disk to be 10 5 — 10 6 s 
(Bath & Pringle 1981). 

There are mainly two types of models for SXTs: the 
disk instability models (DIM) (e.g. Meyer & Meyer- 
Hofmcistcr 1981; Mincshige & Wheeler 1989; Cannizzo 
1992; Cannizzo et al 1995 (CCL)) and the mass transfer 
instability models (MTI) (e.g. Hameury et al 1986: 1987; 
1988; 1990). 

The basic idea for MTI is that the accretion from the 
secondary star is unstable for a range of accretion rates 
roughly from - 10 12 - 10 15 g s" 1 to - 10 16 - 10 17 g s' 1 . 
The secondary star expands under the influence of hard 
X-rays (E > 7 keV) from the primary. The accretion 
rate exceeds the upper limit of the lower stable range 
and jumps to an accretion rate greater than the mini- 
mum of the higher stable range. The subsequent burst 
of mass flow produces the outburst. The disk becomes 
thicker and thicker, finally shielding the region around the 
L\ point. When the X-ray illumination stops the compan- 
ion shrinks and the system returns to its quiescent state. 
MTI models produce the basic characteristics of the light 
curves and the recurrence times. The problem with MTI 
is that, according to the model, the hard X-ray flux with 
E >7 keV must exceed ~ 2.5 x 10 34 (M C /M Q ) 2 erg s" 1 in 
the quiescent state (Hameury et al. 1986) while there is no 
data showing the existence of photons in this energy range 
with such high luminosities before the outbursts (Tanaka 
& Shibazaki 1996). 

DIMs are highly accepted today, although they do not 
produce a completely self consistent picture explaining all 
features of SXTs. Initially DIMs were successful in ex- 
plaining basic characteristics of dwarf novae (Cannizzo 
1994). The model was extended to explain the behav- 
ior of SXTs which are similar to dwarf novae in some 
cases. The disk instability follows a limit cycle mecha- 
nism based on an "S" shaped surface density S versus 
M curve. Upper and lower branches of the "S" corre- 
spond to hot and cool states respectively. The middle 
branch represents an unstable regime. When M > M min 
the disk follows the upper branch switching to the lower 
branch when M < M max . Each radius in the disk has 
its own "S" curve. If the mass transfer satisfies the con- 
dition M m ax(R m ner) < M T < M min {R outer ) then the 
limit cycle mechanism operates through the range between 
Rinner and R ou ter- If this range of radii covers the entire 
disk then the whole disk can jump from one stable branch 
to the other. Accumulation of matter at some radius R 
may cause mass per unit area S(i?) to exceed S TOOX (i?). 
Viscous dissipation increases suddenly. Waves of surface 
density propagate to both smaller and larger radii. As a 
result, the surface density at each radius in the disk ex- 
ceeds the maximum critical value of the lower branch, and 
the disk jumps to the hot branch. At the end of this pro- 
cess the disk finds itself in a high viscosity state. Unlike the 
quiescent state now the viscous time scale becomes small 
and matter flows on to the central object. Because of the 
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matter flow the densities in the outer regions decrease to 
below the critical densities. This causes a cooling front to 
propagate throughout the disk decreasing the local surface 
density at each radius to below the local critical values. 
Propagation of the cooling front means the decay of the 
outburst (e.g., Cannizzo 1992; Lasota et al 1996). 

The main difficulty of this model is that the a values 
needed in order to produce recurrence times of around a 
few ten years are not able to produce the observed am- 
plitude and duration of the outbursts or vice versa (for a 
discussion of the problems see, e.g. Lasota 1996). 

In this work BSXTs will be studied with a different ap- 
proach concentrating on a simple explanation in terms of a 
disk dynamics model for both rise and the decay pattern of 
the light curves until the onset of the first mini outburst. 
We explore the behavior of the disk following a sudden 
release of mass in the outer radius in terms of the disk dy- 
namics model of Lynden-Bell & Pringle (1974) (LP). The 
observations of GS 2000+25 and GS/GRS 1124-68 will be 
used to illustrate the model. 

LP made a study for the evolution of viscous disks in 
general. In particular the solution with no central flux is 
of interest, corresponding to a disk with a black- hole at 
its center. The Green's function solution of LP take ini- 
tial density distributions in the form of delta functions in 
radius. This is an idealization of an initial mass enhance- 
ment in a ring. Convolution of the elementary solution of 
this Green's function model with any initial mass distri- 
bution will give the general solution. Initial mass release 
at a thin ring at the outer radius (or at a specific radius in 
the disk) already corresponds to a solution with a single 
delta function initial mass distribution. Thus the Green's 
function solution of LP is taken here as the physical solu- 
tion for the BSXT outbursts. The motivation leading us 
to this approach is the similarity between the observed X- 
ray light curves of the outbursts and the Green's function 
luminosity calculated by LP. This hypothesis means that 
a sudden dumping of mass in the outer disk leads to the 
observed outburst. The present work applies the model to 
the data from GS 2000+25 and GS/GRS 1124-68. 

The part of the study by LP related to our application 
is summarized in the Appendix. In Section 2 data from 
GS 2000+25 and GS/GRS 1124-68 arc examined with the 
model. The aim is to calculate the photon flux from the 
disk in the observed X-ray band of Ginga (1 — 20 keV) 
as a function of time using the LP model, and compare 
this with observed photon flux data taken with the Ginga 
satellite ( provided by S.Kitamoto, private communica- 
tion). A good lit to the data until the first mini outburst 
is obtained (Figs. 1—4). From this fit the characteristics of 
the outburst, the viscosity parameter a, the total mass re- 
leased and the recurrence times will be obtained for both 
sources. Section 3 summarizes the conclusions, discusses 
the scope of the model and the problems, and tests for 
future applications, including the secondary and tertiary 



maxima commonly seen in the decay phases of the out- 
bursts of BSXTs. 

2. Models For the X-Ray Light Curves of Black- 
hole Soft X-Ray Transients 

LP present a luminosity versus dimcnsionless time, i*, 
curve for the delta function initial condition. This curve is 
similar to the observed light curves of BSXTs. The similar- 
ity is attractive because a mass accumulation at the outer 
disk may also be represented by a delta function initial 
mass distribution. If a disk instability causes a transition 
to a high viscosity state so that the viscous time scale 
becomes short compared to that in the preoutburst state 
then a sudden release of mass would occur. For such a sit- 
uation the Green's function solution with a delta function 
initial mass distribution would approximately describe the 
subsequent evolution of the disk. 

To test this idea we consider a delta function initial 
mass distribution located at a radius R\, T,(Ri,R,0) = 
5m (2-kR)^ 1 5(R — R\) , and use the Green's function in 
Eq.(A18) as the solution of our problem. 5m is the mass 
which is initially injected at the outer radius R\. The 
solution is equivalent to Eq.(A18) except for a factor 5m 
in front. 

We try two viscosity models and compare the results 
to the observations. In the first model we choose the kine- 
matic viscosity v to be constant for simplicity. In the 
second model we use the form of the viscosity parame- 
ter a = ao(zo/R) n where zq is the half thickness of the 
disk. This form is commonly used for the disk instabil- 
ity models. CCL succeeded to produce exponential de- 
cays of the light curves of BSXTs using this form with 
n = 1.5. For the LP model Eq.(AlO) is valid if v is time 
independent, and either independent of R or a power law 
function of R. On the other hand the general prescription 
of the kinematic viscosity is v — (2/3) (P/ p)(a/fl). Since 
(P/p) oc T , v becomes independent of T(R 7 1) only 
when a oc T _1 . This choice corresponds to n = — 2 and 
leads v to be a function of the specific angular momentum 
h as we will discuss in the second model. For both models 
the photon flux is calculated as a function of time, as the 
injected mass at i?i spreads through the disk. The choice 
of a oc T _1 may not be physical. It is implemented 
here in order to apply the LP model. The success of the 
model may indicate that a real disk is effectively repre- 
sented by the LP model with the artificial choices made 
here. Although a temperature and time independent vis- 
cosity is unrealistic, the success of the simple analytical 
LP solution in fitting the light curve motivates us to take 
this prescription as a useful effective representation. We 
shall discuss the viscosity parameter a, the accumulated 
mass 5m and the recurrence times using the two viscosity 
models. 

The models discussed below are summarized in Tables 
(2-4) 
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Table 1. BSXTs with M x > 3M (Ref: Tanaka & Shibazaki 1996 and references therein) 



Name 


BH mass 


Outburst year 




(M) 








(10 35 erg/s) 


(10 15 g/s) 


J0422+32 


> 3.2 


1992 


0.8 


0.2 


0620-00 


> 7.3 


1917,1975 


2 


2 


1124-68 


~ 6 


1991 


7 


7 


J1655-40 


4-5 


1994 


? 


? 


2000+25 


6-13.9 


1988 


3 


3 


2023+33 


8-15.5 


1956,1979 ? 


2 ? 


2 ? 






1989 


0.6 


0.6 



MODEL I : Constant v 

We first investigate a model with v constant in time 
and uniform throughout the disk. We shall obtain limits on 
v by requiring that the (variable) a parameter is less than 
one. If the initial mass distribution is given by Eq.(A.16) 
multiplied by dm the general time dependent solution for 
the couple is Eq.(A18) multiplied by 6m (see Appendix) 



g = 6m 21k 



(x\x) 1 



exp — 



(x\ - xf 



x/x-\ 



(i) 



where I is constant, x = h 1 / 21 = (Sli? 2 ) 1 / 2 ' = 
(GMR) 1 ^ 1 , t, = 2k~H I x\, Fi{z) = e~ z h(z), I,(z) is 
the Bessel function of imaginary argument, and 



K - 2 = 3 {GMf v 



The dimensionless time is 

m ■ 

Substituting 
dh 



t 



in Eq.(A.lO) we have 



- A {GMfvh- 2 



4l 2 K - 2 h 2 -^ . 
Then I = 1/4 for constant v. For < 



(2) 



(3) 



(4) 



(5) 



(g/gi) 
t» 



•C 1 



jx/xi) 



1 



r(i + i) 



(x/xi) 



2U 



(6) 



Substituting I = 1/4, k 2 and /; in Eq.(A.18) we obtain 
the energy dissipation per unit area 



D 



1 



9 



2-kR * OR 
oc i? -3 t* 5 ^ 4 exp 



f -jl + jR/R,) 2 } ] 
\ 2U J ' 



This contains a sharp rise in the exponential factor, and 
decays approximately as a power law at late times. The 
effective temperature of the disk is 



T 



D 



1/4 



, 2 ° 

= c i?- 3 / 4 U 5/16 exp 



-fl- 



8 t» 



(8) 



where a is Stefan-Boltzmann constant and c ~ 3.9 x 
10 7 (6m v i? 2 ) 1 / 4 (cgs). To find the total luminosity radi- 
ated in the disk as a function of time, D must be integrated 
over the surface of the disk, that is, 



L(U) = D 2nR 

JRo 



dr 



For the X-ray luminosity L x to be determined in the 
observational energy band from 1 to 20 keV, the temper- 
ature in Eq.(8) is substituted in the Planck function 



2tt 



h 2 c 2 e e/kT_ 1 



(9) 



and a numerical integration must be performed over the 
X-ray band throughout the disk. Since 

F e = F e [ e, T(U,R)] , 

the luminosity is found for each annular section of the disk 
with thickness dR at a distance R from the center and we 
find the total X-ray luminosity of the disk by integrating. 

The observations report the X-ray photon flux in the 
1—20 keV band rather than the X-ray luminosity. We use 
the photon flux 



N e 



2tt 



h 2 c 2 e e/kT_l 



(10) 



for numerical integration to calculate the total photon flux 
expected from the disk as a function of time, t*. From 
Eq. (8) the temperature of the inner disk is 



(7) To = c R--"^ t 



3/4 ,-5/16 



exp 



-1 

8 



(11) 
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since Rq/R\ <C 1. We designate the dimensionless time 
t* corresponding to the maximum of the light curve by 
i*(max): 



t,=t»(max) 







and we find f„(max)= 0.4. 
T (max) when t* = t* (max) then 

c ~ 1.027 R 3 /4 T (max) . 



If we set T = 



(12) 



Taking the last stable orbit as the inner radius, Rq 
3R g where R g = (2GM/c 2 ), Eq.(8) becomes 



T= 1.027 i?o /4 T (max) R 



3/4 f -5/16 



I 8t. J 



(13) 



T is substituted in Eq.(10) and numerical integra- 
tion is performed over the observed X-ray band (1 — 
20 keV) throughout the disk. 

The total X-ray photon flux n received by the detector 

is 



n 



COS I 

2nd 2 



N 



(14) 



where TV is the total X-ray photon flux radiated by the 
visible surface of the disk, d is the distance of the source 
and i is the inclination angle between the normal of the 
disk and the line of sight of the observer. 



i. GS 2000+25 

The distance of GS 2000+25 is around 2 kpc. The 
black-hole mass is M x = (5.44 ± 0.15) sin~ 3 i M Q 
(Harlaftis, Home & Filippcnko 1996). Taking the in- 
clination angle 47° < i < 75° , the mass range 
is 6.04 < (Mx/M ) < 13.9 (Harlaftis, Home 
& Filippenko, 1996) 

Setting the inner radius at the last stable orbit, Rq = 
3R g , the outer radius of the disk R\ may be chosen ~ 
10 3 -10 4 times the inner radius. Numerical integration for 
photon flux is not sensitive to the outer radius R\ as long 
as Ri ^> Rq, since most of the X-ray radiation comes 
from the inner part of the disk. The maxima of the model 
light curve and the observed light curve can be matched 
by adjusting the single parameter To (max). This must be 
done for different masses tracing the possible M x range. 
The values of T (max) corresponding to different masses 
and different inner radii are given in Table(2). 

Writing the dimensionless time as t* = b t we ob- 
tain the value of the scale factor b that gives the best 
fit (minimum x 2 ) between the model and the observed 
photon count light curves for a grid of mass values in 
the allowed range. We do this after matching the max- 
ima of the theoretical and the observed light curves, since 



there is not enough data during the initial increasing 
phase of the photon flux. The Xmin values have a decreas- 
ing trend with decreasing mass in the mass range of the 
black-hole. Fig.l shows the photon flux versus t* graph 
for M x — 13M Q , with the scale factor b corresponding 
to minimum \ 2 ■ 
For this model 

U = bt = 6 vR~ 2 t. (15) 

Withfo- 5.94 x 10~ 7 for M x = 13 M Q , Eq.(15) gives 

v R\ 2 ~ 8.9 x 10~ 8 s^ 1 . (16) 

In the a disk model (Shakura & Sunyaev 1973) the 
tangential stress can be characterized by the single pa- 
rameter a and the kinematic viscosity, v, may be written 

as 



* (8) 



n 



2 fcT 

3 mpO 



(17) 



Now the question is when and where the viscosity pa- 
rameter a = (3/2) v (mpO/fcT) becomes maximum in 
the disk. Since oc i? _3//2 and kT oc i?~ 3//4 we may con- 
clude that a is always maximum at the inner disk. But 
kT oc 5 ^ 16 implies that the a parameter is increasing 
with time. We follow the luminosity evolution of the disk 
until the time when the inner disk has cooled down to a 
temperature characteristic of the quiescent state. The qui- 
escent state temperatures, or temperature upper limits are 
in the range - 0.1 - 0.3 keV for most BSXTs (Tanaka & 
Shibazaki 1996). We shall use 0.3 keV to derive constraints 
for a. The time corresponding to this temperature is al- 
ready longer than the time of the first mini outburst, after 
which the fit for the decay following the primary outburst 
is no longer relevant. 

Taking kT~ 0.3 keV for the quiescent state inner disk 
temperature and a < 1 we obtain the limits 

v < 1.8 x 10 11 cmV 1 

and using Eq.(16) 

R 1 < 1.4 x 10 9 cm - 95i? - 

i?i may be less than the outermost radius without effect- 
ing the fits as long as the X-ray photon flux coming from 
R > i?i is not significant compared to the R < R\ region. 
R\ denotes the location of mass accumulation and release, 
and is not necessarily identical with the outermost radius. 
In other words we inject 5m at R = R\ in our ideal- 
ized model. We obtain /cTo(max)~ 0.49 keV from the fit 
and c ~ 1.15 x 10 12 (cgs) from Eq.(12). Using the relation 
c = 4.73 x 10 7 {vR^ 2 5m) 1 / ^ (cgs) obtained above for the 
uniform viscosity model, we have 



vR~{ 2 5m ~ 3.5 x 10 17 (cgs) . 



(18) 
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Substituting the value oivR^ 2 from Eq.(16) in Eq.(18) we 
obtain the amount of mass injected initially at the outer 
disk 



6m ~ 3.9 x 10 24 g . 



(19) 



The optically inferred mass transfer rates onto the 
outer disk of BSXTs are around <~ 10 15 — 10 16 g s" 1 during 
the quiescent states. These values are similar to those 
obtained from the long term average, (M x ) , of the 
mass transfer in X-ray outbursts in order of magnitude 
(Tanaka & Shibazaki 1996, and references therein). This 
may show that most of the matter accreting onto the 
disk is accumulated at R\ during quiescence and is used 
as a fuel for the next outburst. This is consistent with 
the observation that M x during quiescence, which reflects 
the mass transfer reaching the inner disk, is much less 
than M pt- Only one outburst was detected for each of 
the sources GS 2000+25 and GS/GRS 1124-68, so we 
use the average accretion rate (M x ) ~ 2 x 10 15 g s _1 
of another typical BSXT, A0620-00, which has a recur- 
rence time around 58 years (White et al 1984), to es- 
timate the recurrence times of our two sources. From 
tx ~ 8m/ (M x ) the recurrence time is found to be around 
62 years for GS 2000+25. 

On the other hand, using a — (3/2) v (mpCl/kT) 
with kT = kT(R out ) ~ 1.67 x 10~ 3 keV where R out ~ 
10 3 x R it is found that a ~ 5.7 x 10~ 3 for the outer 
radius. At R = R±, where mass accumulation occurs, 
a becomes 3.2 x 10~ 2 . These values correspond to an in- 
ner disk temperature of around 0.3 keV. The inner disk 
temperature decreases to 0.3 keV at t <~ 100 days. 

ii. GS/GRS 1124-68 

The distance of GS/GRS 1124-68 is around 3 kpc, the 
black-hole mass is M x ~ 6M Q (McClintock et al 1992). 
fcT (max)= 0.7 keV produces the best fit to the Ginga 
observations of the outburst and decay light curve. For 
the constant v model, analysis similar to the one described 
above for GS 2000+25 gives the results 



vR-{ 



1 x 10 



-7 



5m ~ 3 x 10 24 g 
i rec ~ 49 years 
v < 8.3 x 10 10 cmV 1 
i?i < 9.1 x 10 8 cm ~ 172 R 
a < BSBikT^R- 3 / 2 . 

For GS/GRS 1124-68 we have taken M x ~ 6M Q (McClin- 
tock et al 1992) and from the best fit we have obtained 
b ~ 6.13 x 10~ 7 . When the inner disk temperature re- 
duces to ~ 0.3 keV the a values are ~2x 10~ 2 and 
~ 5.3 x 10~ 3 for R = R x and R = 10 3 i? respectively. 
The best fitting model for the photon flux history is shown 
in Fig. 2 together with the data. The reduced \ 2 is found 
to be 123.7. Although this is a large value, we note that the 



magnitudes of the fluctuations of the observed photon flux 
data about a smooth decay curve are much higher than 
the errors of the observed data points. As shown in Fig. 2 
the model reproduces the average rise and decay behavior 
of the outburst light curve of GS/GRS 1124-68. 



MODEL II: v = v(h) 

In this model, we allow v to vary with radius (equiva- 
lently, with specific angular momentum, h) but not with 
temperature or time. 

i. GS 2000+25 

We adopt the commonly used form of the viscosity 
parameter prescription a = a (z /R) n - Cannizzo et al 
(1995) succeeded to reproduce the exponential type decays 
of the outbursts of BSXTs employing this form of the 
viscosity parameter with n — 1.5. The kinematic viscosity 
is v = (2/3) (P/ p)(a/fl). From the hydrostatic equilibrium 
in the z direction P j p 

4 



fl and zq ~ «s/^ , where 



■ = kT/mp is the sound speed. Then 
2 



v 



(I)" 4n 



(20) 



v will be independent of temperature only for n 
which gives 

2 2 
v = -a (flR 2 ) = -a h . 

The specific angular momentum dependence of v is 
included in the general form of the solution for g given 
by Eq.(l). For this choice of v Eq.(5) becomes 

1 (GM) 2 

2 h~ 



a 



4l 2 



h 2-(l/J) 



(21) 



Then we obtain 
I = 1/3 



a Q (GM) 2 



and 



U = t a (GM) 



1/2. 



R\' 2 



(22) 



(23) 



Using Eqs.(l), (7), (22) & (23) with I = 1/3 we obtain 
the energy dissipation rate 



Doc i, 4,/3 i? 3 exp 



and the temperature 



-[l + (R/R 1 f' 2 



2U 



(24) 



T = ( — 

2(7 



1/4 



= cU^R-V* 



exp 



-{1 + (R/R 1 f' 2 } 
8 £» 



(25) 
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constant 


V 


v = u(h) 


M x 


cos i 


Ro 


fcT (max) 


x 2 ■ 

mm 


kTo (max) 


r ■ 

niin 






(10 6 cm) 


(keV) 




(keV) 




6 


0.26 


5.3 


0.87 


23.5 


0.87 


19.2 


7 


0.26 


6.2 


0.81 


21.9 


0.80 


17.9 


8 


0.47 


7.1 


0.66 


17.8 


0.65 


14.6 


9 


0.53 


7.9 


0.61 


16.3 


0.60 


13.5 


10 


0.58 


8.8 


0.56 


15.1 


0.57 


12.7 


11 


0.61 


9.7 


0.53 


14.2 


0.54 


12.0 


12 


0.64 


10.6 


0.51 


13.4 


0.51 


11.4 


13 


0.66 


11.5 


0.49 


12.8 


0.49 


10.9 



where 



1.1 x 10 11 (5m a i?^ 3/2 ) 1/4 (cgs) . 



(26) 



Now we follow the same procedure as we did for the 
constant v model. The temperature of the inner disk at 
the time corresponding to the maximum of the outburst 
is T (max) . (dTo/dt*) — gives t*(max) = 0.375 and 
T(R,t*) becomes 

T(R,U) ~ T (max) R 3/4 t; 1/3 ir 3/4 

-[l + iR/R.f 2 ] 



x exp 



where R = 3R g = 3(2GM x /c 2 ). 
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(27) 



The corresponding T (max) produces c ~ 1.17 x 
10 12 and, using Eqs.(26), and (29) we find 



6m a R 1 3/2 =i 1.6 x 10 4 g cm" 3 / 2 
5m ~ 3.3 x 10 24 g. 



(30) 



Dividing 5m by M opt ~ (M) ~ 2 x 10 15 g s 1 we obtain 
the recurrence time 

t r <~ 52 years 

The viscosity parameter a, with n = — 2 and z = Vs/Q is 

R 2 ^ 2 

a = cto 2 — (31) 



Table 3. The results obtained from the first model 



with ug = kT/m p and n = (GM x ) 1/2 i?- 3/2 





GS 2000+25 


GS/GRS 1124-68 


R' 2 v- L 


~ 1 x 10' sec 


~ 1 x 10 ' sec 


8m 


~ 3.9 x 10 24 g 


~ 3 x 10 24 g 


tree 


~ 62 years 


49 years 


u(cm 2 s _1 ) 


< 1.8 x 10 11 


< 8.3 x 10 10 


Ri 


< 1.4 x 10 9 cm 


< 9 x 10 s cm 




~ 95i? 


~ 1.7 x 1Q 2 R 


a{R = Ri max ) 


< 3.2 x 10" 2 


< 2.0 x 10" 3 



Once again we perform the numerical integration to 
construct the photon flux of the model disk as a function 
of time. By matching the maxima of the model and the ob- 
served light curve, we obtain different kT (max) values for 
different masses in the range 6 < M x /M < 13.9 . Later 
using i*= bT we obtain b values producing minimum 
X 2 for each mass tested. The results are in Table(2). From 
the best fit we obtain M x = 13M©, R = 1.15 x 10 7 cm, 
fcT (max) = 0.49 keV, and b ~ 4.63 x 10~ 7 . From Eq.(23) 

13 _ D -3/2 ^ 2g ^ 



b ~ 9.34 x 10 13 a R[ 
we obtain 

a R~ 3/2 = 7.29 x 10" 21 cm" 3 / 2 . 



a 



a 



GM x m p 
kT R ' 



(32) 



(29) 



Since kT oc i?~ 3 / 4 the viscosity parameter a cx i? _1//4 
becomes maximum at the inner disk. Then as we have 
found in the previous case a increases with time un- 
til kT ~ 0.3 keV. The inner disk temperature reduces 
to kT ~ 0.3 keV in about 100 days. Eq.(31) with a = 
a max ~ 1 , kT ~ 0.3 keV and R = R = 1.15 x 10 7 cm 
gives a < 1.9 x 10~ 6 . Using Eq.(29) we find 

i?i < 5.3 x 10 9 cm ~ 461i? • 

As stated before we may consider i?i to be the radius 
where the mass accumulates. This Ri value is nearly 5 
times greater than the one obtained in the first model. 

When the inner disk temperature is around 0.3 
keV using Eq.(27) we find that the temperature is ~ 
1.7 x 10~ 3 keV at R ~ 10 3 R and - 3.0 x 10~ 3 keV at 
R = R\ ~ 461 Ro- The corresponding a values are 
0.18 and 0.21 respectively, which are plausible values. 

ii. GS/GRS 1124-68 

For this model kT (max) — 0.7 keV produces the best 
fit to the observed data. Following the steps described 
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100 




Dimensionless Time, t* 



Fig. 1. Photon Flux versus dimensionless time for the constant v model (for GS 2000+25). The crosses are Ginga data points 
obtained from S.Kitamoto. The quoted errors are small compared to the size of the crosses. The circles are extra data points 
we scaled from Fig. 3 in the review paper by Tanaka & Shibazaki (1996). These were not included in fits or the \ 2 calculations. 
They are shown here to indicate the success of the model in encompassing the fast rise and turnover phases of the light curve. 
The dashed line is the LP model fit. 
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"1 1 
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0.1 
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Fig. 2. Photon Flux versus dimensionless time for the constant v model (for GS/GRS 1124-68). 



above for GS 2000+25 we have obtained a similar fit with 
the reduced x 2 = 105.3 (Fig.4). The results are: 

a x R^ 3/2 ~ 8.8 x 10~ 21 cm- 3/2 
5m ~ 2.37 x 10 24 g 
tree ^ 38 years 
a < 1.9 x 10~ 6 
Ri < 3.63 x 10 9 cm ~ 685 R 



a < 2.5 x lQ- 3 (kTRY 



Since M x /i?o = c 2 /6g is constant the upper limit for a is 
the same as that of GS 2000+25 (Eq.32). The recurrence 
time is similar to that of the first model. When the inner 
disk temperature reduces to ~ 0.3 keV the a values are 
~ 0.19 and ~ 0.17 for R = R x and R = 10 3 x R a 
respectively. 
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100 
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"i r~ 

data 
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x 



1 - 



--x-xx- 



0.1 u 1 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 3.5 

Dimensionless Time, t* 



Fig. 3. Photon flux versus dimensionless time for the v = u(h) model (for GS 2000+25). 



From the Roche-lobe geometry the distance from the 
center of the compact object to the inner Lagrangian point 
L\, Rl 1 ~ 2 x 10 11 cm for both sources. When we set 
R\ ~ Rj Jl we find the ctmax *C 1 which is not plausible. 
We may conclude that the radius i?i at which the mass 
accumulation occurs is around or less than 2 x 1CP 2 times 
the disk's outer radius, if the outer radius -R uter ^ s °f ^ nc 
same order as Rl ± . If our inferred 

values of i?i would mean that Rl x is at least fifty times 
greater than i? uter- 



Table 4. The results obtained from the second model 







GS 2000+25 


GS/GRS 1124-68 


5m 




~ 3.3 x 10 M g 


~ 2.4 x lO^ 4 g 


tree 




~ 52 years 


~ 38 years 


Cto 




< 1.9 x 10~ 6 


< 1.9 x 10~ 6 


Ri 




< 5.3 x 10 9 cm 


< 3.6 x 10 9 cm 






~ 461.Ro 


~ 685i? 


a(R = 


Aim..) 


< 0.21 


< 0.19 


a(R = 


10 3 Ro) 


< 0.18 


< 0.17 



3. Discussion and Conclusions 

We have addressed the general characteristics of the out- 
burst light curves of black-hole soft X-ray transients 
(BSXTs), through a study of the typical BSXTs GS 



2000+25 and GS/GRS 1124-68. The Green's function so- 
lution of the accretion disk dynamics model of Lynden- 
Bcll & Pringlc (1974) (LP), with a delta function initial 
mass distribution, produces a rise and decay pattern sim- 
ilar to those of BSXTs. It is most remarkable that this 
simple dynamical model gives both the rise and the sub- 
sequent decay. We have employed this solution adopting 
a delta function initial mass distribution, a mass dm de- 
posited in a ring at radius R\ in the disk, for two different 
time independent viscosity laws. The LP solution restricts 
the viscosity to be a time independent power law function 
of the radial distance from the center of the disk. 

For the first model we chose a constant kinematic 
viscosity v for simplicity. For the second model we 
adopted the commonly used viscosity parameter prescrip- 
tion a = a (z /R) n and chose n = — 2 to remove the 
time varying temperature dependence of v. This leads to 
v = v{h), where h is the specific angular momentum. 

For both models fcTo( max ) ; the inner disk tempera- 
ture at the moment that the light curve for photon count 
rate reaches its maximum, is the free parameter together 
with cos i. We adjusted /cTo(max) to match the maxima 
of the model and the observed photon flux light curves. 
Next, from t* = b t where t* is dimensionless time and t 
is real time, the best fit to the observed photon flux was 
obtained by adjusting the constant b. For GS 2000+25 
this procedure was performed tracing the possible mass 
range, 6M Q < M x < 13.9M , corresponding to the pos- 
sible inclination angle range, 47° < i < 75° (Harlaftis 
at al 1996). For GS/GRS 1124-68 the black-hole mass is 
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Fig. 4. Photon flux versus dimensionless time for the v — u(h) model (for GS/GRS 1124-68). 
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XX 

x X >^> 



0.1 
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_L 



2 3 4 

Dimensionless Time, 



Fig. 5. Photon flux versus dimensionless time for the constant v model for GS 2000+25 (dashed line) together with an expo- 
nential decay curve (dot-dashed line) with a time constant 30 days. The reduced Xmin values are 5 and 12.8 for the exponential 
decay and the model respectively, for fits up to the onset of the first mini outburst (U < 3.5). It is seen that exponential decay 
models for the mini outburst decays will superpose on the continuing decay of the main outburst. By contrast, LP fits to the 
mini outbursts require that each mini outburst "zeroes" the continuing decay of the main outburst, and starts with the same 
initial conditions as the main outburst. 
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M x ~ 6M (McClintock et al 1992). For this source we 
tried four different cos i values (0.25, 0.5, 0.75, 1.0) keep- 
ing M x constant and followed the same steps to obtain 
the best fit. The results are in Tables (2—4) 

For both sources the fits are applied to the data until 
the onset of the first mini outburst. The LP model fits the 
secular decay in the light curve. The large x 2 values reflect 
the excursions about the mean secular decay, which the 
model does not address. Figs. (1—4) show that both mod- 
els reproduce the average characteristics of the outburst 
light curves quite well. The requirement of the model that 
the kinematic viscosity is independent of time and there- 
fore also of the temperature leads to an a parameter which 
increases with decreasing temperature. The usual disk in- 
stability models require the opposite trend of a increasing 
with increasing temperature. Our time independent vis- 
cosity models are artificial and may be considered as an 
effective representation of real viscosity behavior during 
the hot states of the accretion disks of BSXTs. Empirically 
the LP models can explain the rise and the decay following 
the outburst. Our approach here shows that once sudden 
mass release is triggered at some i?i , for example by a disk 
instability, the subsequent light curve can be understood 
as the dynamical evolution of the disk with viscosity. It 
seems plausible that the viscosity mechanism must have 
a temperature dependence at least to start the outbursts. 
The success of the model means that the real disk dynam- 
ics is insensitive to the temperature dependence of the 
viscosity, or that a constant effective viscosity LP model 
is a close approximation to the real luminosity evolution 
with the integrated effects of variable and temperature 
dependent viscosity. To investigate the relation between 
the LP models and disks with variable viscosity is beyond 
the scope of the present paper, and is to be considered in 
subsequent work. 

The mass dumped in the outburst is 5m ~ 10 g 
in all our model fits, for both sources. This is consis- 
tent with accumulation at (M x ) ~ (M x )o620 f° r recur- 
rence times t r <~ 5m/ (M x ) ~ 50 yrs. We expect that 
the mass release is triggered by a disk instability. The 
large amounts of mass release 5m imply surface densities 
S ~ (5m/ R\) ~ 10 4 — 10 5 g cm~ 2 which are far in ex- 
cess of the critical (maximum) surface density values for 
disk instability models, 

S max = 11.4gcm- 2 i?i 05 ^j a- ' 86 

(Shaffer at al 1986). The mass accumulation without 
triggering an instability then requires small a values ~ 
10~ 4 — 10~ 5 in the quiescent disk. At such a the quiescent 
disk could be optically thin. 

There is no complete picture yet which explains all 
properties of BSXTs, in particular the triggering mech- 
anism for the outbursts which drives the accretion disk 
from the quiescent state to the outburst and back to 



the quiescent state. In a recent study by Cannizzo et 
al (1995) exponential decay patterns are produced by the 
disk limit cycle mechanism, but sudden rise of the out- 
bursts and their long recurrence times are not addressed. 

The second (and in two sources, third) outbursts have 
rise and decay patterns which can be fitted with expo- 
nentials with time constants similar to those of the main 
outburst. In addition to this similarity they also pose the 
problem of understanding the burst repetition time (Au- 
gusteijn et al 1993; Chen et al 1993; Mineshige 1994). 
What is the source of the mass producing the second out- 
burst? The data following the second mini outburst does 
not allow a detailed fit with the LP model. Assuming that 
the amplitude of the burst is proportional to the accu- 
mulated mass, 5m, and scaling with the primary outburst 
we find 5m values are — 3 x 10 23 g and ~ 2.5 x 10 23 g 
for the second mini outbursts in GS/GRS 1124-68 and 
GS 2000+25 respectively. Using the time interval between 
the main and the second mini outburst ~ 80 days, we 
find the extra M to supply the 5m of the second mini 
outburst is - 4.5 x 10 16 g s" 1 for GS/GRS 1124-68 and 
- 3.5 x 10 16 g s" 1 for GS 2000+25. This is nearly an or- 
der of magnitude greater than the long term average ac- 
cretion rate, (M). Two possibilities may be considered: 
(1) All the accumulated mass is not released during the 
main outburst. The mass released in the second outburst 
is "leftover mass" and not accumulated between the main 
outburst and the second outburst. (2) The mass flow from 
the companion is enhanced by the X-rays coming from 
the inner disk in the outburst; or both possibilities run 
together. The X-ray heating of the disk may be important 
to trigger the mini outbursts whatever the source of the 
mass creating the second and the third instabilities. 

The time between the main and mini outbursts (~ few 
months)must be a characteristic time scale of the "trig- 
ger" . The response time scale of the upper atmosphere of 
the secondary is too short ( ^ 10 3 s) (Chen et al 1993 
(CLG)) while the viscous time scale is of the order of a 
week. The shielding of the region around the inner La- 
grangian point, L\, has been proposed to explain this time 
interval (CLG). If the geometrical shielding should end 
gradually, the abrupt rise of the first mini outburst at 
a particular time remains unexplained. Two different ap- 
proaches were proposed by Mineshige (1994). An optically 
thick Compton cloud above and below the disk which be- 
comes optically thin in a very short time scale just before 
the first mini outburst may explain the first mini outburst 
by the response of the companion to the X-rays coming 
from the inner disk. Alternatively a second thermal insta- 
bility would be triggered at the outer disk when the strong 
X-ray heating of the disk has raised the temperature and 
decreased the critical density for the trigger sufficiently 
(transient recession of the cooling front). In the former sce- 
nario one still has to explain the secondary mini outburst 
by invoking a different mechanism, and also exhibits sim- 
ilar rise and decay time scales. Similar time interval and 
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characteristics of all the outbursts seem to imply a unique 
mechanism responsible for the triggering of the outbursts. 

It is well known that simple exponential decays after 
the outburst maximum describe the data well. The mini 
outburst decays can also be fitted with exponential decays 
with the same time constants as required by the main out- 
burst(Fig.5). However, the fast rise and the exponential 
decays (FRED) are put together as separate pieces of the 
fit rather than being part of a single dynamical model as 
in the present model. Our LP models fit the main out- 
burst decay quite well. There is an interesting difference 
from FRED models when it comes to trying to fit the 
data incorporating the mini outbursts and their decays. 
In FRED models, the data stream can be fitted well by 
a superposition of FRED models for the main outburst 
and the mini outbursts. For our LP model fits, by con- 
trast, when the fit to the main outburst data is extended 
past the mini outburst (s), it is seen that the model gives 
a count rate that is greater than the observed count rates, 
the deviation starting from the onset of the mini outburst. 
If it is assumed that the ongoing relaxation in the disk af- 
ter the main outburst is stopped, and the onset of the 
mini outburst involves an instability that resets the disk 
to conditions similar to sudden mass release from a local 
accumulation, then the LP model can be applied to the 
mini outburst (s) and their decays. 

Augusteijn et al (1993) drew attention to the simi- 
larity of the main outburst decay, and the mini outburst 
decays. They related this to a feed-back model invoking 
modulation of the mass transfer from the companion by 
irradiation from the outburst. The present approach iden- 
tifies the similarity with repeated conditions of mass ac- 
cumulation and release in the disk itself. The trigger of 
the mass release could be a disk instability. The similarity 
of main outburst and mini outbursts, requires small scale 
repetitions of the main outburst with similar initial con- 
ditions, rather than superposition or convolution with the 
disk state as evolved from the main outburst's decay. The 
presently available data do not allow a detailed fit of the 
mini outbursts with the LP models. 

The behavior of GS 2000+25 and GS/GRS 1124- 
68 are similar to those of A0620-00, and GRO J0422+32. 
The conclusions may be extended to them as well. These 
ideas will be developed and a similar detailed study for 
the other sources will also be attempted in future work. 
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Appendix 

Here we summarize, following Lynden-Bell & Pringle 1974 
(LP), the general arguments of the angular momentum 
and the energy flow in the thin accretion disks. Consider 
a thin annular section of a disk of width dR at radius R 
from the center. The torque or couple acting on the inner 
edge of the annulus is 
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£1 is the angular velocity about the center, v(R) is the 
kinematic viscosity and is the surface density. Vis- 

cous forces are proportional to velocity gradients, and vis- 
cous torques are proportional to angular velocity gradients 
(shear). 

We may write the equation of motion of the annulus 
at radius R in the form 



Wt (hdm) = -M dR - 



(A3) 



where h = flR 2 , dm = Z2irRdR and (D/Dt) = 
(d/dt) + U.V is the total (convective or Lagrangian) 
time derivative following a fluid element. Eq.(A3) states 
that rate of change of angular momentum is given by the 
torque. 

Since D(dm)/Dt = 



is obtained. Lynden-Bell & Pringle (1974) express the de- 
nominator as 



2AvR ^ - 4/ 2 «- 2 h 2 -^ (A.10) 
dR 

where k and 1 are constants. Eq.(AlO) is valid as long as 
v is constant or varies as a function of R, and can be used 
to parametrize the v and h laws. Eq.(A9) reduces to 



dh 2 ~ 4 \i) h df {AA1) 

Resolving Eq.(All) into modes where g oc e~ s * and set- 
ting k 2 — k 2 s 



Dh 



Dt 



dl 
dR 



d 



2itRY,— = --2- = ^ 2nR i Y<v— . 



dR 



dR) 



(A4) 



If the specific angular momentum h = flR 2 is independent 
of time Eq.(A4) becomes 



dg 



d 



dn 



(A5) 



where F is the outward flux of matter through the radius 
R. For a Keplerean disk we can take 



&9 

dh 2 



1 ( k 



h {1/l) - 2 g = (A12) 



and changing the variables x = h 1 / 21 and gi = x l g in 
Eq.(A12) gives a transformation of Bessel's equation 



+ x 



-idgi 



d 2 9l 

dx 2 dx 
with solution given by 



g = e- st (kx) 1 [A(k)Ji(kx) + B(k)J-i(kx)}. (A13) 



n = {GM) 1 ' 2 iT 3/2 



(A6) 



where M is the mass of the star at the center of the disk. 

Energy generated by the matter flowing into the grav- 
itational well is redistributed or convected by gfi, and the 
remainder is dissipated by 



The general solution of Eq.(A.ll) is given by a convo- 
lution of such modes. Then 



exp 



k 2 t 



(kx) 1 [A(k)J t + B(k)J-i] dk 

(AAA) 



D = 



2irR 



dn\ 

dRJ ' 



The continuity equation for the fluid is 



as 

dt 



+ 



dF 



2-kR dR 







(A7) 



(A8) 



where E is the surface density of the disk. Using 
Eqs.(A4&A8) 



where B(k)=0 for the central couple to be zero, g(x = 
0) = 0. A(k) is given by the inverse transform 



/.oo 

A(k)= g(h,0) Ji(kx) {kxf-Ux. (A15) 
Jo 

When there is no central couple as in the case of 
BSXTs any solution may be considered as made up of el- 
ementary solutions whose initial density distributions are 
of the form 



&9_ 

dh 2 



d_ 

dt 



2AvR (dh/dR) 



(A9) 



Y,(Ri,R, 0) = (27ri?i)" 1 5(R - fli) 



(A16) 
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and the corresponding couple is 

g(h u h, 0) = 21 kT 2 xf' 1 5{x - Xl ). (A.17) 

For the special initial condition in Eq.(A17), using 
Eq.(A15) 

A(k) = 21k- 2 x[ k 1 ' 1 Ji{kx x ). 

Substituting A(k) in Eq.(A14) with B(k)=0 

F , (^) (A18 , 

where i* = 2K~ 2 t/x\, Fi(z) — e~ z Ii(z) and I t (z) is the 
Bessel function of imaginary argument. 

For the general initial conditions g(hi,0) the general 
solution is given by 

foe x' {1+l) X 1 

g(h, t) = / g(xf , 0)k- 2 exp - 

Jo t* 

(A19) 

(Lynden-Bell & Pringle, 1974) (LP). Substituting this 
solution in Eq.(A7) and integrating throughout the disk 
one can obtain the total luminosity of the disk for the gen- 
eral initial conditions g(h\, 0). Instead we use the solution 
given by Eq.(A18) which represents a delta function ini- 
tial mass distribution at a radial distance R = i?i , which 
may be considered as the outer radius of the disk. The lu- 
minosity curve obtained by LP with a delta function initial 
mass distribution is very similar to those of BSXTs. 



21k 



X-\ 



exp - 



(xi-xY 



(xi - x) 
2t*x-\ 
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